Energy landscape and rigidity 
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The effects of floppy modes in the thermodynamical properties of a system are studied. From 
thermodynamical arguments, we deduce that floppy modes are not at zero frequency and thus a 
modified Debye model is used to take into account this effect. The model predicts a deviation from 
the Debye law at low temperatures. Then, the connection between the topography of the energy 
landscape, the topology of the phase space and the rigidity of a glass is explored. As a result, we 
relate the number of constraints and floppy modes with the statistics of the landscape. We apply 
these ideas to a simple model for which we provide an approximate expression for the number of 
energy basins as a function of the rigidity. This allows to understand certains features of the glass 
transition, like the jump in the specific heat or the reversible window observed in chalcogenide 
glasses. 

PACS numbers: 64.70.Pf, 64.60.-i, 05.70.-a, 05.65.+b 



I. INTRODUCTION 

The physics of glass formation is a complex multipar- 
ticle problem, and in spite of its importance from the 
fundamental and technological point of view, many im- 
portant questions remain unanswered |lj . As an example 
we can cite the origin of the non-exponential relaxation 
laws or the ability of certain materials to reach the 
glassy state 3] . To tackle these problems there are many 
different approaches 0: phenomenological models like 
the Gibbs-Dimarzio, theoretical theories like the mode 
coupling or the use of extensive computer simulations 
P|. A very interesting question is how the glass transi- 
tion temperature (T g ) depends on chemical composition. 
Chalcogenide glasses (formed with elements from the VI 
column doped with impurities) are very useful for under- 
standing these effects p. As was discovered more than 
2,000 years ago, T g can be raised or lowered by adding 
impurities, and the fragility of the glass can be changed 
from strong to fragile Q- Recently, by using stochastic 
matrices @ , the law that gives the relation between 
T g and the concentration of modifiers [10( has been ob- 
tained, including a constant that appears in the law for 
almost any chalcogenide glass ^lj. Another interesting 
property of glasses is the behavior of their viscosity, which 
is usually referred as fragility .The fragility of a glass is 
also related with the ease of glass formation, in the sense 
that strong glasses are those that do not require a high 
speed of cooling, and fragile glasses are weak glass form- 
ers that require a rapid quench. The ease of glass forma- 
tion can be explained at least in a qualitative way by the 
rigidity theory (RT), introduced by Phillips ^2] an d fur- 
ther refined by Thorpe By considering the covalent 
bonding as a mechanical constraint, the ease of glass for- 
mation is related with the ratio between available degrees 
of freedom and the number of constraints. If the number 
of constraints is lower than the degrees of freedom, there 
are zero frequency vibrational modes called floppy |14| . 
The resulting network is under-constrained. A transition 
occurs when the lattice becomes rigid. Glasses that are 
rigid due to a certain chemical composition are easier to 



form, and many features of this transition have been ex- 
perimentally observed 6] 16] . Even for simple systems 
like hard-disks 01 an d colloids ^^|, it seems that rigid- 
ity plays an important role . For more complex systems 
like proteins, rigidity has been used as a very powerful 
tool to understand folding and long-time scale motions 
0- 

A very puzzling fact of RT that has not been explored 
is the following: according to the idea of looking at rigid- 
ity as a vectorial percolation problem, at the rigidity 
threshold the entropy is high [21|, due to strong fluc- 
tuations as happens in any phase transition. One even 
can define a free energy and specific heat as a function 
of the flexibility of the system, that has a singularity at 
the transition [22j. However, the experimental data from 
modulated scanning calorimetry in chalcogenide glasses 
shows the opposite: at the rigidity transition the configu- 
rational entropy is less and there is a window of reversibil- 
ity 6] 23]. Specially, it has been observed that protein 
foldin g is reversible because it occurs at the rigidity tran- 
sition I13|. and this seems to be a crucial property for life 
to exist [l9j . 

Mainly, the problem resides in the fact that although 
RT provides a framework to understand many features of 
a system, its use in a quantitative way has not been fully 
developed to provide a link with the thermodynamics of 
the system 20] . In a previous paper we approached this 
problem by using a phenomenological free energy to ac- 
count for many thermodynamical properties of the glass 
transition |2-'1|, and then we made extensive computer 
simulations with associative fluids to show that many 
concepts of the RT work in a "thermodynamical envi- 
ronment" [25||26j. However, the connection with ther- 
modynamics is still not mature, since there is no general 
way of introducing thermodynamics in the RT. 

In a different context, the energy landscape is a formal- 
ism that has been very useful for describing the molecular 
scale events that happen during the glass transition [27| • 
The landscape is a multidimensional surface generated by 
the system potential energy as a function of the molecu- 
lar coordinates Q . In an TV body system the landscape is 
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thus determined by the potential energy function, given 
by $(r 1; r^v) where Yi comprise position, orientation 
and vibration coordinates. For the simplest case of par- 
ticle possessing no internal degrees of freedom, the land- 
scape is a (3N + 1) object. The topography of this land- 
scape is fundamental for the thermodynamics of the sys- 
tem. At high temperatures the system does not feel the 
summits and valleys of $(ri, ...,rjv) because the kinetic 
energy contribution dominates. However, as the temper- 
ature is lowered the system is unable to surmount the 
highest energy barriers and therefore is forced to sample 
deep minima. When this happen, the kinetic of relax- 
ation changes from exponential to stretched exponential 
P|. An important observation is that, according to sta- 
tistical mechanics, the entropy of the system depends on 
the accessible volume in the phase space. However, inside 
a local minimum of the potential energy, it can happen 
that if there are no paths that connect to other minima, 
the system cannot sample that part of the phase space. 
In such a case, ergodicity is broken and the system is 
no longer in thermal equilibrium. Such a glass will have 
a residual entropy |2§|. In this article, we show that 
rigidity can be related with the statistics of the energy 
landscape, since the number of floppy modes is related 
to the number of different configurations of the system 
with nearly equal minimal energies, and thus provides 
an estimation for the number of minima energy basins of 
the landscape. But floppy modes also provides channels 
in phase space that increase the entropy, which in part 
explains the paradox of the window of reversiblity. To 
show these connections, we will concentrate in the effects 
of rigidity on the shape of the energy landscape. 

The layout of this work is the following: in section II we 
discuss a simple way to introduce thermodynamics into 
RT, however, as we will see, the straightforward manner 
of doing this do not agree with the experimental results. 
Thus we propose that the effects of floppy modes are only 
important at low temperatures or during glass transition. 
In section III the connection with the energy landscape 
is made and a simple model is worked out. Finally, in 
section IV we give the conclusions. 



II. RIGIDITY AND THERMODYNAMICS 

In this section we explore some simple thermodynam- 
ical consequences of the RT. As explained before, the 
rigidity ideas of Phillips [3l| and Thorpe were used 
in order to understand the ease of glass formation. In this 
theory, the ability for making a glass is optimized when 
the number of freedom degrees, in this case 3N, where 
N is the number of particles, is equal to the number of 
mechanical constraints (N c ) that are given by the bond 
length and angles between bonds. 

The number (3 A'' — N c )/ 3N gives the fraction of cyclic 
variables of the Hamiltonian, i.e., when one of these 
variables is changed, the energy of the system does not 
change, as for example happens with the center of mass 



coordinate. This fraction also corresponds to the frac- 
tion of vibrational modes with zero frequency (/), called 
floppy modes, with respect to the total number of vibra- 
tional modes. The counting of floppy modes in a mean- 
field, known as Maxwell counting, goes as follows |29| : 
since each of the r bonds in a site of coordination r is 
shared by two sites, there are r/2 constraints due to dis- 
tance fixing between neighbors. If the angles are also 
rigid, in 3-D there are (2r — 3) constraints, to give, 

f _ 3N-N c _, ^ [r/2 + (2r-3)]x r _ 5, 
J ~ 3A ~ ^ 3 6 {T) 

r 

where the last term corresponds to the angular con- 
straints, x r is the fraction of particles with coordination 
r, and (r) is the average coordination number, defined 
as, 

(r) = rx r 

r 

A rigidity transition occurs when / = and the system 
pass from a floppy network to rigid one. If / is a negative 
number, i.e., if there are more constraints than degrees of 
freedom, the lattice is overconstrained and the important 
number is how many stressed bonds are present. In 3D, 
the rigidity transition leads to the critical value (r c ) = 2.4 
if all angular constraints are considered. In real systems, 
the Maxwell counting breaks near the rigidity transition, 
and the number of flop py modes are obtained form the 
pebble game algorithm |14j . 

What are the simple thermodynamical effects of floppy 
modes? To answer this question, first we use the most 
simple model for atomic vibrations in the harmonic ap- 
proximation, where the interatomic potentials are re- 
placed by springs. The corresponding Hamiltonian is, 

3N p2 3N{l-f) 
J=l 3=1 

where Qj and Pj is the j— normal mode coordinate in 
phase space, and uij is the corresponding eigenfrequcncy 
of each normal mode. Observe that floppy modes have 
zero frequency; they do not contribute to the elastic en- 
ergy. Using simple statistical mechanics, we can obtain 
the partition function in the canonical ensemble at the 
classical limit (high temperatures compared with the De- 
bye temperature), 




where T is the temperature and k the Boltzmann con- 
stant. The free energy of the system is now given by, 

3AfcT (2-nmkT\ kT 3N ^ f \ ( 2nkT\ 
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From this last expression, the corresponding specific heat 
(CV) is, 

In this simple approach, the prediction is that Cy is given 
by the Dulong-Petit law, minus a term that depends on 
the number of floppy modes. The reason is clear: floppy 
modes do not store energy since they are cyclic variables 
of the Hamiltonian, as for example happens with the cen- 
ter of mass. However, a careful examination of the ex- 
perimental data shows that for chalcogenide glasses [(|, 
like compounds of As — Ge — Se, and ferroelectric ma- 
terials 15], Cv does not depends on /. Instead, they 
follow the Dulong-Petit law. From this simple thermo- 
dynamical argument, one is lead to propose that floppy 
modes do not have a perfect zero frequency, i.e., in real 
glasses they are shifted by residual forces, like the Van 
der Waals interaction. This argument is confirmed by 
neutron scattering experiments, where it has been shown 
that floppy modes in As — Ge — Se are blue-shifted 0] , 
forming a peak at around hmeV. Thus, at high temper- 
atures, all the 3iV oscillators are excited. Only at low 
temperatures we suggest that the effects of floppy modes 
are important since all floppy modes are frozen nearly 
at the same temperature. The corresponding tempera- 
ture (6/) where these modes are frozen, can be estimated 
from the energy required to excite modes of 5meV, that 
gives 6/ ~ 60°^. 

This behavior at low temperatures, where a quantum 
treatment is needed, can be modelled by using a simple 
density of states p(ui) that takes into account the floppy 
peak in the spectrum. First we use a Debye type of den- 
sity of states, normalized to 3jV(1 — /).Then we add the 
contribution from the floppy modes, with a delta func- 
tion centered around a characteristic peak at uiq. The 
corresponding density of states is, 



9JV(1-/) , ,2 
3 ^ 



3Nfd(uj — luq), if lo < luo 
if uj > ljd 



where ujd is the Debye cut-off frequency. By using the 
Bose-Einstein distribution for the number of phonons in 
equilibrium at a certain temperature, we get that the 
specific heat is, 

C v = (1 - f)3NhD(x ) + f3Nk 



(e* - l) 2 

where x — Qfi/T, xo = 8d/T, and &d = Tvud is 
the Debye temperature. D(xq) is the well known De- 
bye function. At high temperatures, the model predicts 
the Dulong-Petit law as expected, while at low T, the 
following behavior is obtained, 



Cv ~ (1 - f)3Nk 



4tt 4 ( T 
67 



contribution is determined from the fraction of floppy 
modes for a given composition of the glass. The present 
model suggest that experiments at low temperatures per- 
formed on chalcogenide glasses will provide caracteristic 
features of rigidity. 



which is a Debye law of the type T 3 , but with a contri- 
bution that is in the form of the Einstein model. Each 



III. ENERGY LANDSCAPE AND RIGIDITY 



In the last section we discussed that floppy modes have 
effects mainly at low T . In spite of this, an examina- 
tion of the experimental results shows that the number 
of floppy modes is also important for the thermodynami- 
cal properties at the glass transition 0- For example, the 
magnitude in the jump of Cp, usually denoted by ACp, 
the jump in the thermal expansion, the energy for activa- 
tion of viscosity, the fragility and the entropy of a liquid 
melt depends on /. Moreover, very recently the group of 
Boolchand discovered the window of reversibility in the 
heat flow, associated with a phase of zero internal stress 
in the lattice [2j| . Angell has pointed out the qualitative 
relationship between energy landscape and fragility dur- 
ing glass transition [2]|. However, still is no clear how to 
relate these features with the statistics of the landscape. 
Here we will show that the number of floppy modes pro- 
vides a useful parameter to represent the roughness of 
the landscape. This roughness is evident when the glass 
is melted explaining why floppy modes are important 
during glass transition, since they are collective motions 
that provide pathways across the phase-space and energy 
landscape. 

As a first and tentative step, we start again by suppos- 
ing that floppy modes are at zero frequency. Around any 
given inherent structure, the potential has a minimum 
and thus can be expanded in a Taylor series, which turns 
out to be the expression of an harmonic potential. From 
the Hamiltonian presented in eq.Q is clear that in a in- 
herent structure, each floppy mode provides a channel in 
the landscape since the energy does not depend upon a 
change in a floppy coordinate. A very simple example 
is shown in fig. 1, which shows the bottom of the land- 
scape for a system with two normal modes. In the first 
system, (Fig. la) / = 0, but the other has / = 1 (Fig. 
lb) since one of the springs constants was set to zero 
(of course, by excluding the center of mass coordinate). 
In a more general way, for a given inherent structure, 
the number of channels is clearly given by /. Each chan- 
nel increases the available phase space allowed to visit. 
The entropy due to floppy modes is easy to calculate. In 
the microcanonical ensemble, the number of accessible 
states (£l(E, V, N)) for a system with a volume V is pro- 
portional to the area defined by the surface of constant 
energy, E = H(P 1 , ...Pn, Qi, ...,Qn). Since floppy modes 
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are cyclic variables of the Hamiltonian, we can write, 



:!A ! 



E=H(P 1 



3JV(1-/) 
,J=1 fc=l 



and using the Boltzmann relation S 
get, 



fclnf2(£,V,iV) we 



5 = In 



{2 7 rm)—E 3N ( 1 



h™ [(3JV(1 - //2)) 




(2) 



The entropy provided by the channels in the landscape 
is simply given by the last term, S c = fNk In V. At first 
glance, it seems that this result agrees with the exper- 
imental observations, because during glass transition, it 
has been observed that floppy glasses have a great en- 
tropy and as a result, they have a more fragile behavior 
as deduced from the Adams- Gibbs relation Q. However, 
a more detailed analysis shows that if we suppose an en- 
tropy of the type given by eq. J5J, the specific heat does 
not follow the Dulong-Petit law. This is due to the de- 
pendence of S upon E 3N ( 1 ~f/ 2 \ which is just a result 
of the independence of H with respect to floppy modes. 
As discussed in the previous section, this leads to the 
conclusion that floppy modes are not strictly at zero fre- 
quency. The blue-shift of the floppy modes means that 
the channels in phase space are not flat: there is a small 
curvature in the direction of the floppy variable. This 
effect has the property that it restores the Dulong-Petit 
law and provides directions in phase space where the sys- 
tem can relax without big changes in energy. 

In a floppy glass there is a hierarchy in the strength 
of the forces. The forces that restores the Dulong-Petit 
law are the weakest. Then it is natural to assume that 
the anharmonic contributions of these residual forces are 
also small. Under this assumption, the extra entropy due 
to these modes is S ~ fNklnV which is only activated 
when the glass traverses the glass transition. Further- 
more, we can speculate that these channels are in fact 
the ones that explains the fragility and ease of glass for- 
mation since is clear that is much more difficult to trap 
the system in a local minima of the landscape when many 
channels are present. 

However, there are two important facts to consider in 
all the previous statements: first the number of floppy 
modes is a function of the energy. In fact, when the glass 
becomes fluid, most of the constrictions upon the bond 
lengths and angles are relaxed and / is raised. For the 
extreme case of no-bonding between atoms, the system 
behaves without constraints and all the modes are floppy 
/ = 1. Notice that an ideal gas is a perfect "floppy sys- 
tem" . An improvement to eq. (2J is to make / a function 
of E, then the number of floppy modes is 3Nf(E).ln such 
a case, the jump in the specific heat will also depend on 
/, as observed in the experiments. The function f(E) 
is zero when E 3> kT g and has a value determined by 



the average coordination number below the glass transi- 
tion, i.e., f(E) = 2 — | (r) . The shape of this function 
v^can'be^stimated using a procedure that we will describe 

The second consideration is that the number of floppy 
modes affects the number of minima energy valleys (usu- 
ally called inherent structures) that are available when 
the system has a certain energy. This effect is explained 
in figure 2, where a system of bars and hinges is consid- 
ered. In the example of figure 2, there are no angular 
/iVfcln^ ces ' Each bar provides a restriction to the system. 
There are three squares. In one of the squares there 
is a diagonal bar. As a result, this square can not be 
deformed, since the distance between all the hinges are 
fixed. The other two squares are flexible as indicated 
by the arrows. Each of these flexible squares can be de- 
formed independently, and the system has 2 floppy modes 
(again, without counting the center of mass translation 
and rotations around it). Now we move the diagonal 
bar to the second square and the system has the same 
number of floppy modes, but the structure is different, 
and the same thing happen if we put the diagonal in the 
first square. In the landscape formalism, each of these 
configurations is in a different "inherent structure" and 
corresponds to a basin with the same energy. This part 
of the entropy has been studied extensively in the context 
of rigidity transitions pl| . However, as we will see next, 
there is a competition between the channel and configu- 
rational entropies. 

To see how these concepts are applied in a particu- 
lar case, let us consider the following two dimensional 
model that contains all of the previous features that we 
discussed. Consider a system of N disks interacting with 
a central force where no angular forces are considered. 
Each disk has a hard core potential and an attractive 
part which has a range determined by the parameter A. 
If a is the diameter of the disks, r is the distance between 
the centers of two disks, the potential is written as, 



3C 



V{r) = 



-Vi 







if r < a 

if a < r < Act 

if r > Act 



The nature of the fluid and solid phase of this system 
has been studied in a previous work |18(. Here we only 
study the rigidity. Within this model, a bond is formed 
when the distance between two disks is between a and 
Act. Each bond has an energy — Vi, and the energy of 
the system is just proportional to the number of bonds. 
This number is proportional to the average coordination 
number divided by 2 since each bond is shared by two 
sites. Then, the amount of energy (E) of the system is 
given by, 

E = -ViA^ * -2V 1 N(l - /) (3) 

where it was used that for the mean field approxima- 
tion in two dimensions / ~ (27V - (N < r > /2))/2N. 
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From the last equation, is clear that a gas is obtained 
when the system is 100% flexible (/ = 1) and the state 
of maximal packing (the hexagonal lattice with maximal 



coordination . 



6) is overconstrained (there are N/2 



redundant bonds in the mean field approximation). No- 
tice that / is a function of E. 

As said previously, there is an entropy provided by 
floppy modes channels (Si) and to the different configu- 
rations of floppy modes (S2). According with our previ- 
ous assumptions, the first contribution is Si ~ fNklnA 
where A is the area of the system. This is only valid in 
the flexible phase, i.e., before the freezing of the system 
since at that point it has been suggested that there is a 
rigidity transition pfij . After freezing, this contribution 
is zero (Si = 0). At high temperatures, the system is a 
fluid and the entropy is just the same as the one obtained 
from the available phase space without any interaction. 
A more realistic assumption although still very rough is 
to use that Si ~ fNk \n(A — b) where b is proportional 
to the area occupied by the disks [3i| b « Nn(\a) 2 /2. 

The other contribution to the entropy comes from the 
number of ways in which a configuration with a given (r) 
can be made. Although this number is difficult to calcu- 
late, one can suppose a cell model of the fluid, and then 
just consider the number of ways in which absent bonds 
can be deleted from the lattice with maximal packing. 
This number of configurations (Q(f,N)) is, 



tl(f,N) 



(r)N 
2 



(r)N 
2 



(4) 



where (r) is a function of /. The corresponding config- 
urational entropy is S2 = khiQ(f, N). A natural way 
to compute this entropy is to define an order parameter 
m(f) as, 

m(/)= (r--(r»-(r) B V f l (5) 

In terms of this parameter, and using Stirling's approx- 
imation, the total entropy for / > now reads, 



Si + S 2 
Nk 



ln2 + fhi(A-b) 
(m(f) + 1) 



1 ft , m\ (1 - m(f)) 
ln(l + m(f)) ln(l 



For an overconstrained lattice, the expression for the en- 
tropy is just S2 The expression for / > contains the 
effects that were discussed previously, i.e., the linear de- 
pendence of the entropy upon /, and the contribution 
from different structures with the same energy. In figure 
3 we show a plot of the total entropy and the correspond- 
ing contributions for a given A — b. It is interesting to 
note that S2 tends to grow as we diminish the number 
of floppy modes, since the number of configurations with 
the same energy grows. Notice that 5*2 does not have a 
maximum exactly when / = due to the mean field ap- 
proximations; the maximum is shifted to the right. From 



eq. |0, this occurs near (r) = 4, .i.e., near the two di- 
mensional rigidity transition. This fact seems to contra- 
dict that in the rigidity transition, the experimental non- 
reversible heat flow is a minimum, which means that the 
configurational entropy is a minimum. One can expect 
that in the rigidity transition, a lot of fluctuations will be 
observed, while in the experiments it seems that the con- 
trary is true [3 ■ However, the present results shows that 
floppy modes have two competing effects, one is the en- 
tropy due to the different configurations, but the other is 
the shape of each basin, since around inherent structure, 
floppy modes form channels that increase the entropy. 
Thus, as is shown in figure 3, when the system pass from 
flexible to rigid, the number of configurations raises, but 
the number of channels diminish . Experimental results 
suggests that this last effect is more important, since the 
configurational entropy of a melt with a floppy glass for- 
mer is higher as the number of floppy modes is increased 

00 

Finally, the free energy for / > can be written as, 
W- = -2Vi(l-f)^-T(Si + S 2 ) (6) 

To compare with the energy landscape formalism, we 
use that the partition function is the sum of partition 
functions at inherent structures 1341, 



Z(T) = Z ha (T) f G(E)e~ E/kT dE 
Jo 



(7) 



where G(E) is the number of energy basins with energy 
E, and Z^ a (T) is the partition function for a system of 
harmonic oscillators 34]. Since G(E) is always a growing 
function, and e~ E ^ kT is always decreasing, the integral 
of eq.(JZJ) can be replaced by the value at the maximum 
E, 



G(E)e~ E / kT dE a G(E)e- E ' kT 



The corresponding free energy F is, 



^ = -\nZ(T) = -\nZ ha (T) - InG(E) + A (8) 



As usual, the free energy is just the contribution from 
rn 6fe 1 Vibrations inside the basin, the entropic component 
due to the existence of different basins, and an energetic 
component which reflects the average depth of the land- 
scape at a certain T . Now we turn our attention in how 
G(E) is affected by the floppy modes. Comparing eq.@ 
and eq.©, we get an estimation for G(E), 



G(E) « expiV 



-m(E) \n(A - b) - ( m ( g ) + ^ ln(1 + m (E)) 



where m(E) is obtained from eq.@ and eq.JSJ, 
mTO-l-VCBJ-^ + jgy) 



(9) 
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If the channel term is the most important, G(E) can be 
approximated by, 

G(E) « exp [N \m(E)\ (In A - (b/A))] 

which has the same general shape of that proposed by 
Stillinger [35|. The factor \m(E) (In A - (b/A)) in the 
exponential can be identified with the landscape com- 
plexity j3(|. 



IV. CONCLUSIONS 

In this article, we have explored the effects of floppy 
modes into the thermodynamics of glasses. In particular, 
we shown that a blue-shift of floppy modes can be pre- 
dicted using simple thermodynamical arguments. This 
leads to the formulation of a simple model, which sug- 
gests effects of floppy modes at low temperatures. Dur- 



ing glass transition, floppy modes also play a role. Thus 
we explored how flexibility and rigidity determine the 
energy landscape. We found two competing effects that 
contribute to the entropy in the liquid melt; one contribu- 
tion is given by channels and the other is the existence of 
different energy basins. By considering a simple example, 
we showed how to estimate both contributions, and we 
discussed the effects in the window of reversibility. The 
results of this article seem to confirm the Phillip's idea 
that glass forming tendency is enhanced at the rigidity 
transition [37| . since although there is an increase in the 
entropy due to the different energy basins, the pathways 
provided by floppy modes are absent and the system is 
easier to trap in a certain minimum. 
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Figure 1. Bottom of the landscape for a system a) 
with no floppy modes and a potential energy given by 
4>(x, y) — x 2 + y 2 b) with one floppy mode obtained by 
removing a "spring" , 4>(x, y) — x 2 .A channel is generated 
in the y direction. 

Figure 2. A system of bars and hinges with three dif- 
ferent configurations. The squares with the diagonal bars 
are rigid, while the others are flexible. The corresponding 
floppy modes are shown with arrows. 

Figure 3. Contributions to the total entropy (crosses) 



in units of Nk. The dotted line is the contribution from 
channels (Si) with the aribitrary value A — b = 4. The 
solid line is the contribution from different configurations 
(S2). 
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